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Abstract 



A new model for the finite one-dimensional harmonic oscillator is proposed based upon the 
algebra u(2)q,. This algebra is a deformation of the Lie algebra u(2) extended by a parity 
i-S^ . operator, with deformation parameter a. A class of irreducible unitary representations of u(2)q, 

is constructed. In the finite oscillator model, the (discrete) spectrum of the position operator 



(~| . is determined, and the position wavefunctions are shown to be dual Hahn polynomials. Plots 



of these discrete wavefunctions display interesting properties, similar to those of the parabose 
oscillator. We show indeed that in the limit, when the dimension of the representations goes 
to infinity, the discrete wavefunctions tend to the continuous wavefunctions of the parabose 
oscillator. 

m 
> , 

5 ■ 1 Introduction 

V/^ ■ Finite oscillator models obey the same dynamics as the classical and quantum oscillators, but the 
operators corresponding to position, momentum and Hamiltonian are elements of some algebra 
different from the traditional oscillator Lie algebra. The interest in finite oscillator models comes 
primarily from optical image processing [l]. In this context, most attention has been paid to a 
finite oscillator model based on the Lie algebra su(2) (or so (3)), generalized to so (4) in the case of 
a finite two-dimensional oscillator [1^3]. 
^ ' In the one-dimensional case, the setting is as follows. There are three (essentially self-adjoint) 

■ operators: a position operator q, its corresponding momentum operator p and a (pseudo-) Hamilto- 
nian H which is the generator of time evolution. These operators should satisfy the Hamiltonian-Lie 
equations (or the compatibility of Hamilton's equations with the Heisenberg equations): 

[H,q] = -ip, [H,p]=iq, (1) 

in units with mass and frequency both equal to 1, and h = 1. Contrary to the canonical case, the 
commutator [q,p] = i is not required. Wigner considered such a system already in 1950 [1]. He 
required the extra condition 

H=^(p^ + f). (2) 

In that case, one is dealing with the Wigner quantum oscillator (or parabose oscillator) OE]. In 
fact, the algebraic structure equivalent with ([1]) and ([2]) is just the Lie superalgebra osp(l|2) (see 
Appendix). It turns out that in this case one is still dealing with an oscillator with an infinite 
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energy spectrum, the only difference being that this spectrum is shifted compared to the canonical 
case. This shift is determined by the osp(l|2) representation parameter a (a is real and positive), 
the spectrum being n + a (n = 0, 1, 2, . . .). 

If one wants a finite oscillator model, the relation ^ should be dropped and one is left with ([T]) 
only. On top of ([T|) and the self-adjointness, it is then common to require the following conditions [2]: 

• all operators q, p, H belong to some (Lie) algebra (or superalgebra) A] 

• the spectrum of H in (unitary) representations of A is equidistant. 

The first requirement is often taken to be stronger, in the sense that all commutator brackets 
between operators close into a Lie algebra (or its universal enveloping algebra) [2]; it is also often 
taken to be weaker, in the sense that the remaining commutator [q,p\ is just required to be a 
"function" of [7j. Also the second requirement is sometimes dropped [7]. 

The case with A = 5u(2) (or its enveloping algebra) has been treated extensively in a number 
of papers [THS]. In that case, the relevant representations are the common su(2) representations 
labelled by an integer of half-integer j. Up to a constant, the Hamiltonian H is the diagonal 
5u(2) operator with a linear spectrum re + ^ (n = 0, 1, ... , 2j). Obviously also q and p have a 
finite spectrum. In this model it is given by + [2j. More interestingly, the 

position wavefunctions have been constructed, and are given by Krawtchouk functions (normalized 
symmetric Krawtchouk polynomials). These discrete wavefunctions have interesting properties, 
and their shape is reminiscent of those of the canonical oscillator [2]. It was indeed shown that 
under the limit j — )• oo the discrete wavefunctions coincide with the canonical wavefunctions in 
terms of Her mite polynomials [2l[8]. 

In the present paper, we construct a one-parameter deformation of the enveloping algebra of 
u(2), with a real deformation parameter a > —1. This algebra, denoted by u(2)q, is chosen as 
the algebra A to which the operators q, p and H should belong, along with the relations ([1]). For 
the value a = — ^, u(2)q reduces to the undeformed u(2) (or su(2)). For half-integer j values (not 
for integer j- values), the common su(2) or u(2) representations can be extended to representations 
of u{2)a- An interesting aspect of these representations, is that the operators q and p have again 
a simple spectrum. Furthermore, the position wavefunctions can be computed explicitly. These 
wavefunctions turn out to be dual Hahn functions (normalized dual Hahn polynomials). The dual 
Hahn polynomials, usually characterized by two parameters (a,/3), appear in the wavefunctions 
with parameters (q, a -|- 1) or (a -|- 1, a), depending on the parity of the wavefunction. We present 
some plots of these discrete wavefunctions, which have interesting properties similar to those of the 
5u(2) finite oscillator model (to which they reduce when a = — ^). To understand these properties, 
the limit j — )• oo is determined. Quite surprisingly, under this limit the discrete wavefunctions 
become the parabose oscillator wavefunctions. So, the finite oscillator model presented here could 
be interpreted as a finite parabose oscillator model, which reduces to the su(2) finite oscillator 
model when a = — ^. 

The structure of the paper is as follows: in Section 2 we construct the algebra u(2)q, and its 
representations. In Section 3, it is shown that u{2)a gives rise to a new finite oscillator model. 
The main result here is the determination of the spectrum of the position operator q, and the 
computation of its eigenvectors. In Section 4, we present the position wavefunctions of the u(2)q, 
oscillator model. We discuss some of their properties, and give the limit relation to parabose 
wavefunctions. Finally, we have included an Appendix where some aspects of the parabose oscillator 
(or Wigner quantum oscillator) are summarized. 
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2 The algebra u(2)q^ and its representations 

The Lie algebra 5u(2) [HulOj is usually defined by its basis elements Jq, J+, J_ with commutators 
[Jo, J±] = ±J± and [J+, J_] = 2Jq. The non-trivial unitary representations of su(2), corresponding 
to the star relations Jq = Jq, J± = Ji^ , are labelled [9l[T0] by a positive integer or half- integer j. 
These representations have dimension 2j + 1, and the action of the su(2) operators on a set of basis 
vectors \ j, m) (with m = —j, + +j) is given by 



Jo\j^^) = ^ \j,rn), J±\j,m) = =F ± m 1) \j,m± 1). 

The Lie algebra su(2) can be extended to u(2) by adding an extra operator C which commutes with 
all basis elements Jo, J+, J-; the action of C in the above representations is diagonal: C\j,m) = 
(2j -|- 1) |j, m). It is also possible to extend u(2) further by a parity operator P, whose action in 
these representations is given by P\j,m) = (—1)-?+™ \ j,m). On the algebraic level, this means that 
we extend the universal enveloping algebra of u(2) by an operator P that commutes with C and 
Jq, that anticommutes with J_|_ and J_, and for which P^ = 1. 

The addition of these two operators does not yet give much extra structure as far as represen- 
tations are concerned. But this extended algebra can be deformed by a parameter a, leading us to 
the definition of u{2)a- 

Definition 1 Let a be a parameter. The algebra u(2)q, is a unital algebra with basis elements Jq, 
J+, J_, C and P subject to the following relations: 

• C commutes with all basis elements, 

• P is a parity operator satisfying P^ = 1 and 

[P, Jo] = PJo - JoP = 0, {P, J±} = PJ± + J±P = 0. (3) 

• The su(2) relations are deformed as follows: 

[Jo,J±] = ±J±, (4) 
[J+, J_] = 2 Jo - (2a + ifP - {2a + l)CP. (5) 

Relation ([5]) looks at this moment fairly arbitrary. It would also be of little interest, if the repre- 
sentations of u(2) could not be extended to representations of u(2)q; for the current choice of ([5]), 
this happens to be the case. Note that for a = — ^, ([5]) reduces to the usual su(2) commutator 
relation; so for a = the deformation is trivial. 

The extension of the common (2j -|- l)-dimensional representations of u(2) to representations of 
u{2)a is possible provided j is a half-integer (and not when j is an integer). This is presented in 
the following proposition. 

Proposition 2 Let j be a half-integer (i.e. 2j is odd), and consider the space Vj with basis vectors 
|j, — j), \j, + . . ., |j, j). Assume that a > — 1. Then the following action turns Vj into an 
irreducible representation space ofu{2)a. 

C\j,m) = {2j + l)\j,m), (6) 

P|j,m) = (-iy+™|j,m), (7) 

MJ^'m) = m\j,m), (8) 



j \/(i - m){j +m + l) \j, m + 1), if j + m is odd; 

\ ^/ {j — m + 2a + + m + 2a + 2) | j, m + 1), if j + m is even. 



\ + m){j - m + l) \ j, m - 1), if j + m is even; 

J-\j,m) = < (10) 

\\/{j + m + 2a + l){j — m + 2a + 2) \j,m — 1), if j + m is odd. 
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Note that, for a general parameter a, the condition 2j odd is necessary in order to have J+|j, j) = 
and J-\j,—j) = 0. The condition a > —1 ensures that the factors under the square roots are 
positive (a = — 1 would lead to reducibility). To actually prove the proposition it is sufficient to 
verify that all defining relations from Definition 1 are valid when acting on an arbitrary vector | j, m), 
and this is a straightforward calculation. Irreducibility follows e.g. from the fact that {J^)^\j, — j) 
is nonzero and proportional to \ j, —j + k) for A; = 1, 2, ... , 2j, and similarly {J-)^\j,j) is nonzero. 

Note also that the representation given in this proposition is unitary under the star conditions 



3 A one-dimensional oscillator model based on u{2)a 



The conditions for finite oscillator models have been described in Section 1. In the context of the 
algebra u(2)q, one can make the following choice for the position, momentum and Hamiltonian 
operator: 



1 



( J+ + J- 



P ■ 



h = Jq + Ic. 



It is easy to verify that ([T]) is satisfied. Furthermore, in the representation space Vj, H\j,m) 
{m + j + m), so the spectrum of H is indeed linear and given by 



1 

^+2 



(n = 0,l,...,2i). 



(12) 



An interesting question is the determination of the eigenvalues of the operators q (position 
eigenvalues) and p (momentum eigenvalues), and the corresponding eigenvectors. We shall perform 
this for the position operator (for p, the analysis is similar). 

From the actions ([^- (|1U|) . one finds 

2q\j, m) = \/ {j + m){j - m + 1) \ j, m - 1) + {j - m + 2a + + m + 2a + 2) \ j, m + 1) 
if j + m is even, and 

2q\j, m) = Y^(j + m + 2a + 1) ( j — m + 2a + 2) | j, m — 1) + ^/ {j — m){j + m + 1) \ j, m + 1) 



if j + m is odd. This means that in the (ordered) basis B oi Vj, 



B = {\j,-j),\j,-j + 1) 
the operator 2q takes the matrix form 



2q 



with 

M = f V(fe + l)(2i-fc), if is odd; 

^ [ ^J{k + 2a + 2){2j - A; + 2a + 1), if k is even. 

For this matrix, the eigenvalues are known explicitly |lHll2j. 



(13) 
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Ml 








Ml 























M2j-i 



(14) 



(15) 
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Proposition 3 The 2j + 1 eigenvalues q of the position operator q in the representation Vj are 
given by 

11 1 , , 

- a- j - -,-a - j + -,. . . ,-a - l;a + l,a + 2,... ,a + j + -. (16) 

It will be appropriate to label these q-eigenvalues as q^, where k = —j, — j + 1, . . . , +j, so 

= ±(a + + -), k = 
In the 5u(2) oscillator model [2], the spectrum of q in the representation Vj with 2j odd is given by 

-j, -J + 1, • • • , ^, • • • , j + 

So in the u(2)q, oscillator model, there is a shift from the origin by a + ^ (so for positive g-eigenvalues 
a shift by a + ^; and for negative g-eigenvalues a shift by —{a + ^)). Note that for a > —\ this 
shift is away from the origin; for — 1 < a < — ^ this shift is towards the origin. 

In a recent paper [12], it was shown that also the eigenvectors of a tridiagonal matrix of the 
form (|14|) can be constructed explicitly. The expressions of these eigenvectors involve Hahn poly- 
nomials, so let us first recall some notation. Hahn polynomials Qn{x; a, N) [TSKH] of degree n 
(n = 0, 1, . . . , N) in the variable x, with parameters a > —1 and /3 > — 1 are defined by [l3l[T4] : 

n ^ » AT\ T? ( -n,n + a + P + l,-x \ 



in terms of the generalized hyper geometric series 3F2 of unit argument |15lll6j . Hahn polynomials 
satisfy a (discrete) orthogonality relation [13]: 



N 



J2 a, /3, N)Qi{x; a, /3, N)Qn{x; a, /?, N) = h{n; a, /3, A^) (18) 



x=0 



where 



w{x;a,f^,N)=(^ ^ j j {x = 0,1, . . . , N); 

w o (n + a + /3 + l)jv+i(/3 + l)„n! 

h[n; a, p, N) — 



{2n + a + l3 + l){a + l)„(iV - n + 1)„7V! ' 



We have used here the common notation for Pochhammer symbols |15yi6j {a)k = a{a + 1) • • • (a + 
k — 1) for k = 1,2, .. . and (a)o = 1- As it; is the weight function and h{n;a, (3, N) the "squared 
norm", orthonormal Hahn functions Q are determined by: 



Qn{x; Q, /3, N) ^ : . (19) 

y/h{n;a,l3,N) 

Dual Hahn polynomials are almost the same as Hahn polynomials, but with the role of x and n 
in (fT7|) interchanged. So for x G {0, 1, . . . , A^}, the right hand side of p!7|) can be seen as the dual 
Hahn polynomial of degree x in the variable A(n) = n(n + a + /3 + 1) [13] . 
The next result follows now from |12t Proposition 2]: 

Proposition 4 The orthonormal eigenvector of the position operator q in Vj for the eigenvalue qt, 
denoted by \ j,qk), is given in terms of the basis B by 

j 

\j,qk)= Yl Uj+m,j+k\j,m). (20) 
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Herein, U = {Urs)o<r,s<2j is a (2j + 1) x (2j + 1) matrix with elements 

U2r,j-s-^- = U2r,j+s+l = ^^^-(^^ a,a + l,j- (21) 
?72,+i,,_,_i =-C/2,+i,,+,+ i =-^a(r;a + l,a,i-^), (22) 

where r, s G {0, 1, . . . , j — The functions Q are normalized Hahn polynomials (|19p . 

Note that U is an orthogonal matrix, UU'^ = U'^U = I, hence the q eigenvectors are orthonormal: 

{j,Qk\j,qi) = Ski- 

4 u{2)a oscillator wavefunctions and their properties 

In general, the wavefunctions are the overlaps between the normalized eigenstates of the position 
operator and the eigenstates of the Hamiltonian. So the wavefunctions of the u(2)q, finite os- 
cillator are the overlaps between the g-eigenvectors and the //-eigenvectors (or equivalently, the 
Jo-eigenvectors |j, m)). They are denoted by ^^•+'^(9), where m = —j, —j + 1, . . . and where q 
assumes one of the discrete values qk {k = —j, + +j). Concretely, following the notation 

of the previous section: 

'^f+milk) = U, m\j, qk) = Uj+rn,j+k- (23) 

Let us examine this function in more detail. Since U has a different form for even and odd 
indices, (j2ip and (j22p . we shall also make this distinction here. For j + m even, j + m = 2n, this 
is by ([2T|l an even function of the position variable q; for positive g-values one has 

^2n(.Qk) = — ^Qfe_i(ra;a,a + 1,J - -), k = -,-,... ,j; 
or more explicitly: 



^(a). ^ (-1)" / w{n;a,a + l,j - ^) f -qk + a + l,qk + a + I, -n \ 

^2n Wfc) = —^\ T7 : —r-. — TT 3i^2 1 ; 1 • (24) 

V2 \l h{qk-a-l;a,a + l,j - ^) \ a + l,-j + ^ J 

So one can interpret this as a dual Hahn polynomial of degree n in the position variable q^ (or 
rather in the variable X{k — ^)). For j + m odd, j + m = 2n + I, this is by ([22|) an odd function of 
the variable q; for positive g-values we have 

(a) (~1)" ~ 113 

$2n+i(%) = ^/|-Qfe_i(^;a + l,a,J - k = -,-,..., j; 

or explicitly: 

f . (-1)" / w{n;a + l,a,j - i) f -qk + a + l,qk + a + 1,-n \ 

again a dual Hahn polynomial of degree n in the position variable. Because of the appearance of 
these polynomials as wavefunctions of the u(2)q finite oscillator, we shall refer to this model as the 
Hahn oscillator. 

It is interesting to study some plots of these discrete wavefunctions, for some values of a. Let 
us choose a fixed value of j, say j = and plot some of the wavefunctions ^l^\q) for various 
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values of a. Since a = — ^ is a special case (where u(2)a reduces to u(2)), there are three cases 
to be considered: — 1 < a < — i, a = — | and a > — |. In Figured] we take a = — |, a = —0.7 

and a = 1 respectively. We also plot in each case the ground state (^Q^\q), some low energy states 

$["''(g) and ^^2\l)^ 9''^'^ the highest energy state. 

For a = these plots are familiar. In that case, the current finite oscillator model coincides 
with the model of Atakishiyev et al [HE], based upon the su(2) algebra. It is known that in such 
a case, the wavefunctions ^^^^(9) are in fact Krawtchouk functions. This is indeed a special 
case of our wavefunctions, expressed as dual Hahn functions. This follows from the fact that when 
a = — i the dual Hahn polynomials, which are 3F2 series appearing in reduce to 2F1 

series according to 



(26) 



/-g + l/2,g + l/2,-n \ Qj) /-2n,-j-q_ 

^^^i 1/2, -,+1/2 'V"^"'^ (^py^'^^l -2j ' 

f-q + l/2,q + l/2,-n \ (-1)" (an+i) f -'2n - 1, -j - q \ 

'"^"K 3/2, -,+1/2 'V-"^(-:p)^'^^i -2j ''J- ^'"^ 

These reductions have been given in [1^ and can be obtained, e.g., from [H (48)]. The 2F1 series in 
the right hand side correspond to symmetric Krawtchouk polynomials (i.e. Krawtchouk polynomials 
withp = 1/2 [13]). 

Another reason why the wavefunctions <I>i ^^'^\q) can be considered as finite oscillator wave- 
functions is because in the limit j — t- 00 they yield the ordinary oscillator wavefunctions. To present 
this limit, one should also pass from a discrete position variable g to a continuous position variable 
X. This can be done by putting q = j^^'^x, and then taking the limit. As a consequence of this factor 
j^/^, due to the discrete orthogonality relation for the functions $n ^^'^\q)^ which should become a 
continuous orthogonality relation in the limit case, one should consider the limit of j 
Then one obtains: 

lim = -^^^H4x)e-^y^ (28) 



2"/2Vn!7rV4 

where Hn{x) are the common Hermite polynomials |13^ ll7j. For the current case with a = — ^, 
the discrete wavefunctions are Krawtchouk functions and this limit has been computed in [T]. 
Expression (j28p can also be deduced from our limit for general values of a. 

Now the main question is: what is the j — )• cxd limit of (^^\q) for general a {a > —1)? To 
compute this, let us consider the case n even and n odd separately. For even values, one should 
consider expression (f24]l . and compute 

limji/4$H(//2^). (29) 

Since q is positive in (|24p . x is also positive; the total wavefunction is even so it is clear how to 
extend it to negative x-values. The limit of the 3F2 function in (f24l) is quite easy: 



1- T? ( -3 ' x + a + l,3^/^x + a + l,-n \ f -n 2\ r(»)f 2n fnn\ 

lim 3i^2 ^ , . ^ 1 ; 1 = iFi • X = -— — >{x ), (30) 

j^oo y Q + l,-J + i J \a + l J {a + l)n 

where is a Laguerre polynomial |13pi7j. So it remains to determine 



V2 ^ h{j^/^x-a-l;a,a + l,j -^)' ^ ^ 
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In order to compute this, one should replace all Pochhammer symbols in the expression of h by 
Gamma-functions, i.e. (a)/. = T{a + k)/T{a). Then it is a matter of combining appropriate factors 
in the quotient of (j3ip . and using Stirling's approximation for the Gamma-function [17j. This leads 
to: 

lim j^^^w{n;a,a + l,j ^ {a + l)n ^2a+i^-x\ (32) 

j^oo 2h{j^/'^x - a - I; a, a + 1, j - ^) n!r(a-|- 1) 

Combining all these results, and extending it trivially to negative x-values, one finds: 



lim //^cI>£)(//2,) = (-W |.r+V2,-.V24«)(,2), (33) 

j-5>oo y T{a -I- n -I- 1) 



For odd values of n, the limit of the 3F2 function in (|25p yields: 

( -j^/'^x + a + l,j^/'^x + a + l,-n \ ^ / 

(34) 

Performing a similar computation as before (and extending it to negative x-values as an odd 
function of x), the final result is: 



3^00 y L [a + n + 2) 

Note that for a = —1/2, one finds indeed (f28]l . 

The functions in (jSSp and (j35p are familiar: they are in fact the wavefunctions '^^^^\x) of the 
parabose oscillator with parameter a = a -|- 1 > (see the appendix). So we have: 

lim j'^'^tHj'^'^) = ^t^'Hx). (36) 

In this sense, the current model can be interpreted as a finite one-dimensional parabose oscillator 
model. This also explains the shape of the discrete wavefunctions plotted in Figure 1. For — 1 < 
a < the shape typically reproduces the continuous wavefunctions of the parabose oscillator 
with 0<a = a + l<^: see the plots for a = —0.7 in Figure 1 and those for a = 0.3 in Figure 2. 
For a > — ^, the shape of the wavefunctions is different than for a < Now the shape is similar 
to those of the parabose oscillator with a > | : compare the plots for q = 1 in Figure 1 with those 
for a = 2 in Figure 2. 



5 Discussion 



The quantum harmonic oscillator in the canonical non-relativistic case has a very simple and 
well known solution, with equidistant energy spectrum and stationary states described in terms 
of Hermite polynomials. It is used as a model in many applications. However, it also has some 
restrictions. In particular, it has an infinite spectrum (which is sometimes not realistic), and its 
wavefunctions satisfy a continuous orthogonality relation with infinite support. Therefore, it is not 
directly applicable to describe models where only a finite number of eigenmodes can exist, such as 
in optical image processing. 

In this context [1] one realized that the quantum oscillator equations of motion are consis- 
tent with commutators other than [q,p] = i. Fractional Fourier transforms [18, Ch. 10] to signal 
analysis on a finite number of discrete sensors or data points led to physical models realizing a 
one-dimensional finite oscillator |19fl21j . 
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Finite oscillator models are specific examples of quantum systems with a finite Hilbert space. 
Such finite quantum systems have been studied from a general point of view, see the review pa- 
per [22] and references therein. A wide variety of applications of finite quantum systems is known, 
of which quantum optics and quantum computing are the most popular. Since our contribution 
is devoted to a model of the one-dimensional finite oscillator, primarily applications in quantum 
optics are worth considering. Having this in mind, it is necessary to turn to a two-dimensional 
finite oscillator model and study wavefunctions on a finite (square) Cartesian grid [2]. In the sim- 
plest picture, this follows just from a direct product of two one-dimensional oscillators, and the 
wavefunctions in an xy-plane could be described, in an obvious notation as in [2], as 

These Cartesian mode wavefunctions could be depicted as in [2, Figure 3], with however two 
main differences. First of all, the "density distribution" depends on the a-parameters, as in the 
one-dimensional plots of Figure 1. Secondly, the "sensor points" of the grid are not uniformly 
distributed, but according to the position spectrum given in Proposition 3. We hope that these 
extra facilities open the way to more sophisticated techniques for two-dimensional signals on square 
screens, in the spirit of [23]. This analysis, however, falls outside the scope of the present paper. 

The most interesting (mathematical) aspects of finite oscillators were given in [IHS] . The model 
of the finite oscillator considered there is based on the Lie algebra 5u(2) and its representations 
labelled by an integer or half- integer j. The authors have shown that the wavefunctions of the su(2) 
finite oscillator are given by Krawtchouk polynomials, and that in the limit j — )• oo these discrete 
wavefunctions tend to the canonical oscillator wavefunctions in terms of Hermite polynomials. The 
position operator has a finite equidistant spectrum. These properties are summarized in the bottom 
part of Table 1. 



finite u(2)q oscillator 
(Hahn oscillator) 
algebra u(2)q, 
discrete position spectrum with a gap 
dual Hahn polynomial 



parabose oscillator 
(Wigner quantum oscillator) 
Lie superalgebra osp(l|2) 
position spectrum: M or M \ {0} 
Laguerre polynomial 



I (a 



J, (a 



or a 



finite 5u(2) oscillator 
(Krawtchouk oscillator) 
Lie algebra su(2) 
discrete equidistant position spectrum 
Krawtchouk polynomial 



canonical oscillator 
oscillator algebra 
position spectrum: R 
Hermite polynomial 



Table 1: Summary of some properties of the four oscillator models appearing in this paper. For each 
of them, we give the dynamical algebra, the spectrum of the position operator, and the polynomials 
appearing in the wavefunctions. The arrows indicate how to go from one model to another. 

In the present paper we have been able to extend the 5u(2) finite oscillator model by introducing 
an extra parameter a > —1. The dynamical algebra, u{2)a is a deformation of the Lie algebra u(2) 
extended by a parity operator. The representations of u(2)q, are those of u(2) characterized by a 
half-integer j, but deformed by the parameter a. The spectrum of the position operator has been 
determined: it is again finite and equidistant in steps of one unit, except that there is a gap of size 
2a -|- 2 in the middle of the spectrum. The wavefunctions of the u(2)q, finite oscillator have been 
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determined, and turn out to be dual Hahn polynomials involving this parameter a. When q = — |, 
the su(2) finite oscillator model is recovered, depicted in the left part of Table 1. 

We have also investigated the j — )• oo limit of the u(2)q finite oscillator model. It is quite 
remarkable that under this limit, the discrete wavefunctions in terms of dual Hahn polynomials tend 
to the continuous wavefunctions of the parabose oscillator |24II25| (in terms of Laguerre polynomials, 
sometimes referred to as generalized Hermite polynomials). This correspondence is given in the top 
part of Table 1. Note that this parabose oscillator, often referred to as the one-dimensional Wigner 
quantum oscillator, is itself an extension of the canonical oscillator. The parabose oscillator is 
described [26] in terms of the Lie superalgebra osp(l|2), with (infinite-dimensional) representations 
labelled by a positive number a (in our correspondence o = a + 1). In the representation with 
a = ^, the parabose oscillator coincides with the canonical oscillator. 

Our introduction and study of the Hahn oscillator has not only extended the by now well known 
5u(2) finite oscillator model, it has also unified two approaches that extend the canonical oscillator. 
In this unification, given schematically in Table 1, both the parabose oscillator model and the 5u(2) 
finite oscillator model appear as a limiting case or a special case. 

Our approach is based on a deformation of the Lie algebra u(2) and its representations. This 
deformation is from a very different nature than the common g-deformation suq{2) in the context 
of quantum algebras: in that case, the common dynamical algebras become quantum algebras, and 
the wavefunctions are deformed into the corresponding g-functions. Note that suq{2) has also been 
considered to build a model for a finite oscillator [27l|28]. In that case, the position operator has 
a discrete anharmonic spectrum, and the wavefunctions are given in terms of dual g-Krawtchouk 
polynomials |28j . 

Note, finally, the relation with a recent paper on the shifted harmonic approximation [29]. In 
that paper, the authors define in a sense discrete wavefunctions as Clebsch-Gordan coefficients of 
SU{2) or SU{1, 1). They study some plots of these discrete functions, and also show that in the 
limit these functions approach the harmonic oscillator wavefunctions. Since SU{2) or SU{1,1) 
Clebsch-Gordan coefficients can also be expressed as 3i<2(l) series, i.e. as Hahn or dual Hahn 
polynomials [30^31j. the link with the current paper is clear, giving in a way a more fundamental 
reason why Clebsch-Gordan coefficients can be considered as discrete wavefunctions. 

We hope to investigate in a more systematical manner how the algebra u(2) can be deformed, 
as in ([5]), in a way that is still consistent with operator actions in representations. Such a general 
approach might lead to a unification at a higher level, in terms of polynomials at a higher level in 
the Askey-scheme. 



Appendix 

Let us summarize some aspects of the parabose oscillator, or the one-dimensional Wigner quantum 
oscillator (see also |25j for an overview). Consider the Hamiltonian for a one-dimensional harmonic 
oscillator: 

where p and q denote respectively the momentum and position operator of the system. Wigner [3] 
already noted that there are other solutions besides the canonical one if one only requires the 
compatibility between the Hamilton and the Heisenberg equations (dropping the canonical com- 
mutation relation [q,p] = i). These compatibility conditions are: 

[H^p\=iq, [H,q] = -ip. (38) 
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So, one has to find operators f) and acting in some Hilbert space, such that the compatibiHty 
conditions (p8|) hold, with H given by ([37|) . Also, in this Hilbert space p and q have to be self- 
adjoint. 

The solutions to (I38p can be found by introducing two new operators 6^ and 6^ (the parabose 
creation and annihilation operators): 

b^ = ^{q^ip), (39) 

or equivalently 

q = l=(h^ + b-), p=^ib+-b-). 

It is then easily checked that 

H = l{b-,b+}, (40) 



and that the compatibility conditions (|38p are equivalent with 

[{6-,6+},6±]=±26±. (41) 

The relations (j4ip are in fact the defining relations of one pair of parabose operators b^ |32] . 
Moreover, from the self-adjointness of the position and momentum operators it follows that 

(6±)t = 6^. (42) 

It is known that the Lie superalgebra generated by two odd elements b^ subject to the restric- 
tion (|^T]) is the Lie superalgebra osp(l|2) [26]. So, the solutions to the problem are given by the 
star representations of the Lie superalgebra osp(l|2). These are known, and are characterized by a 
positive real number a and a vacuum vector |0), such that 

6-|0) = 0, {b-,b+}\0) = 2a\0). 

The representation space can then be shown to be the Hilbert space ^^(Z+) with orthonormal basis 
vectors \n) (n S Z+) and with the following actions: 

b+\2n) = v^2(n + a)|2n + l), 6-|2n) = \/2^ |2n - 1), 

b+\2n + 1) = x/2(n + 1) \2n + 2), b~\2n + 1) = ^/2{n + a) \2n), 

from which it immediately follows that 

{b-,b+}\n) =2{n + a)\n). (44) 
The energy spectrum follows from (j44p and (j40p : 

H\n) = (n + a) 

so one has an equidistant energy spectrum with ground level given by a. This, and the explicit 
action of the commutator action [p,q] on basis vectors \n), confirms that only a = 1/2 yields the 
canonical solution. 

To construct the wave functions, one can work as follows. Using q = (6+ + 6~)/\/2, one finds 
now from (j43l) 

q\2n) = -^{V2^\2n - 1) + y^2{n + a)\2n + 1)), 



q\2n + 1) = ^(V2(n + a)\2n) + ^2{n + l)|2n + 2)). (45) 
v2 
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This means that q is (or extends to) an unbounded Jacobi operator on This Jacobi operator 

corresponds to the generahzed Hermite polynomials, and the spectrum of g is M |33j . Consider the 
formal eigenvectors v{x) of q, for the eigenvalue x, and write these again as 



The equation qv{x) = xv{x) leads, using (|45p . to a set of recurrence relations for the coefficients 
^n^\x). Solving these recurrence relations explicitly, taking into account the normalization condi- 
tion, one finds, in terms of Laguerre polynomials: 



Since these coefficients have an interpretation as the position wavefunctions of the Wigner oscillator, 
we are done. Alternatively, one can work in the position representation, where the operator q is still 
represented by "multiplication by x", and the operator p has a realization as —i^ plus an extra 
term (depending on a) [341 Chapter 23]. Using this realization the time- independent Schrodinger 
equation can be solved, also yielding the expressions (|47l) [24]. 

To compare with the discrete case considered earlier, let us also plot in Figure [2] the wavefunc- 
tions '^^n\x) for n = 0, 1, 2 and for a = 0.3, a = \ and a = 2 (corresponding to a = —0.7, a = — | 
and a = 1 respectively). Note that the shape of the discrete parabose oscillator wavefunctions are 
similar to those of the continuous parabose oscillator, a fact explained by the limit relation. 
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Figure 1: Plots of the discrete wavefunetions ^l^\q) in the representation with j = 65/2. The 
four top figures are for a = —1/2, the middle figures for a = —0.7, and the four bottom figures for 
a = 1. In each case, we plot the wavefunetions for ra = 0, 1, 2, 65. 
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Figure 2: Plots of the parabose oscillator wavefunctions (x). The three figures on the left are 
for a = 1/2 and correspond to the canonical case; the figures in the middle are for a = 0.3, and 
three figures on the right are for a = 2. In each case, we plot the wavefunctions for n = 0, 1, 2. 
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